Pressure Induced Valence Transitions in /-Electron Systems 
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A review is given of pressure induced valence transitions in /-electron systems calculated with the 
self-interaction corrected local spin density (SIC-LSD) approximation. These calculations show that 
the SIC-LSD is able to describe valence changes as a function of pressure or chemical composition. 
An important finding is the dual character of the /-electrons as either localized or band-like. A 
finite temperature generalisation is presented and applied to the study of the p-T phase diagram of 
the a — > 7 phase transition in Ce. 

PACS numbers: 



I. INTRODUCTION 

The knowledge of valence of rare earth or actinidc 
ions is important for the understanding of the solid state 
physical properties of /-electron systems. It determines 
for example the equilibrium volume of /-electron com- 
pounds. This is extremely well illustrated in work of 
JayaramaniS. In Fig. Q] we show the lattice constants 
for the rare earth sulphides, selenides and tellurides (af- 
ter Ref. @). The abrupt expansion of the lattice for 
SmS, SmSe, SmTe, EuS, EuSe, EuTe, TmSe, TmTc, 
YbS, YbSe and YBTe is associated with a change in 
valence from trivalent rare earth to divalent rare earth. 
Hence this figure shows a direct correlation between lat- 
tice constants and valence. Whilst the determination of 
the valence is important for the static properties of /- 
electron systems, it becomes less useful for understand- 
ing the dynamical properties such as heavy fermion and 
Kondo screening of local spin magnetic moments. How- 
ever it was shown 3 " that the construction of an alloy anal- 
ogy for the valences can describe the finite temperature 
behaviour, in particular the pressure versus temperature 
phase diagram of Ce. Therefore dynamical fluctuations 
between valences could possibly describe the finite tem- 
perature behaviour of such systems as heavy fermions. 

The /-electron systems are a rich hunting ground 
for structural and valence transitions as a function of 
pressureA^ Advances in experimental techniques, such 
as progress in high pressure cell technology, make pres- 
sure experiments feasible over an ever-increasing range, 
possibly up to 1000 kbar. Furthermore the constant im- 
provements in brightness of the synchrotron sources al- 
low for the study of high pressure phases of /-electron 
materials in greater detail. 

In this paper we will review a methodology for ab 
initio calculations of valence and hence valence transi- 
tions as a function of pressure. The methodology is 
based on the self-interaction corrected local spin density 
(SIC-LSD) approximation. SIC-LSD corrects for a spu- 
rious self-interaction of an electron with itself. This self- 



interaction is in most circumstances small and hence the 
LSD is sufficiently accurate for most applications. How- 
ever for localized states this correction is substantial and 
the SIC-LSD approach needs to be applied. The SIC- 
LSD method differentiates between localized and itiner- 
ant electrons and leads to an orbital dependent potential. 
Hence, self-interaction corrections capture a dual picture 
of coexisting localized and band-like /-electronsi^iS As 
a function of volume or pressure rearrangements in the 
number of localized and band-like /-electrons take place 
and a study of this forms one of the topics of this paper. 

The paper is organized as follows. In Section II the 
basic aspects of the SIC-LSD method are outlined, con- 
centrating on its implementation in terms of bands and 
k-points. In Section III we comment on the so-called lo- 
cal SIC implementation and its generalization to alloys 
and finite temperatures, allowing to study both static va- 
lence and spin fluctuations. Section IV contains the cal- 
culations of valence of elemental rare earths, their mono- 
sulphides and mono- nitrides. Also the results for the 
valences for Ce-, Sm-, Eu- and Yb-mono-pnictides and 
mono-sulphides are presented in that Section. Section 
V discusses the valences of the actinides and their com- 
pounds. The results for the a — > 7 phase transition in 
Ce at finite temperatures is the subject of Section VI. 
Section VII presents the summary and conclusion. 



II. THE SIC-LSD TOTAL ENERGY METHOD 
AND VALENCE 



The total energy functional of the local spin density 
(LSD) approximation to density functional theory pro- 
vides good accuracy in describing conventional solids 
with weakly correlated electrons However, for all 
rare earths and for materials containing actinidc elements 
(beyond Np) the electron correlations are significant, and 
the self-interaction correction needs to be included to ob- 
tain a similar accurate description of the localized nature 
of the 4/ or 5/ electrons. 
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The ensuing SIC-LSD total energy functional! 3 , is de- 
rived from the LSD as: 



^SIC-LSD 
j^LSD 



pLSD 



T + U + E xc + V ext +E so , 
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where a labels the occupied electron states and <5f 
is the self-interaction correction for state a. As usual, 
E LSD can be decomposed into a kinetic energy, T, a 
Hartree energy, U , the interaction energy with the atomic 
ions, Vexti an d the exchange and correlation energy, 
E xc i^ Here in addition we include the spin-orbit cou- 
pling term E so . The spin-orbit energy, for each occupied 
state a is: 



(5) 



The A.E S i C term expresses the spurious self-interaction 
of the LSD energy, which is subtracted from the LSD 
energy in equation (1). The self- interaction is calculated 
for each occupied state a, and is given by the sum of the 
Hartree interaction and the exchange-correlation energy 
for the charge density of that state: 
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E T 
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For itinerant states, 5^ IC vanishes identically, while for 
localized (atomic-like) states the self-interaction may be 
appreciable. Thus, the self-interaction correction con- 
stitutes a negative energy contribution gained by an /- 
electron when localizing, which competes with the band 
formation energy gained by the /-electron if allowed to 
delocalize and hybridize with the available conduction 
states. The volume dependence of 8 a is much weaker 
than the volume dependence of the band formation en- 
ergy of rare earth 4/ or actinide 5/ electrons, hence the 
overbinding of the LSD approximation for narrow / band 
states is reduced when localization is allowed. The SIC- 
LSD energy functional in Eq. (1) appears to be a func- 
tional of all the one-electron orbitals, but can in fact be 
viewed as a functional of the total (spin) density alone, 
as discussed in Ref. 0. 

One major advantage of the SIC-LSD energy func- 
tional is that it allows for different valence scenarios to be 
explored. By assuming atomic configurations with differ- 
ent total numbers of localized states, self-consistent min- 
imization of the total energy leads to different local min- 
ima of the same functional, E SIC ~ LSD in Eq. (I), and 
hence their total energies may be compared. The config- 
uration with the lowest energy defines the ground state 
configuration. Note, that if no localized states are as- 
sumed, E SIC ~ LSD coincides with the conventional LSD 
functional, i.e., the Kohn-Sham minimum of the E L 
functional is also a local minimum of i^sic-lsd _ Among 



the actinide elements, this was also found to be the global 
minimum for U and Np, but for the later actinide ele- 
ments, and also in all rare earths, minima exist with a 
finite number of localized states^ The reason is that the 
respective / orbitals are sufficiently confined in space to 
benefit appreciably from the self-interaction correction. 

Another advantage of the SIC-LSD scheme is the possi- 
bility to localize /-states of different character. In partic- 
ular the various crystal field eigenstates, either magnetic 
or paramagnetic. 

The SIC-LSD still considers the electronic structure of 
the solid to be built from individual one-electron states, 
but offers an alternative description to the Bloch picture, 
namely in terms of periodic arrays of localized atom- 
centered states (i.e., the Heitler-London picture in terms 
of Wannier orbitals). Nevertheless, there still exist states 
which will never benefit from the SIC. These states retain 
their itinerant character of the Bloch form, and move in 
the effective LSD potential. This is the case for the non-/ 
conduction electron states in the rare earth and actinide 
metals. 

Results of two implementations of the SIC-LSD scheme 
will be presented in this paper. The first approach ap- 
plies the SIC to the localized bands, whilst in the second 
approach the SIC is applied to the phase shifts of the 
localized states. Both approaches can be considered to 
be a reasonable ansatz, for the solid state, to the imple- 
mentation of the SIC for atoms. The latter approach is 
local in nature - as expressed by the use of single site 
phase shifts- whilst the former can take into account a 
more extended behaviour of the wave function. How- 
ever the band-based approach involves repeated back and 
forth transformations of the wavefunctions from recipro- 
cal space to real space localised functions. In real space 
the band-dependent SIC potential is evaluated, in recip- 
rocal space the bandstructure problem is solved. This 
is a time consuming aspect of this method. Further de- 
tails of this implementation can be found in Ref. E For 
a variet y of applica tions to d and / electron solids, see 
Refs. I7ll8ll9l20l2ll and references therein. The results 
in Sections IV and V are based on applying the SIC to 
the localized bands, whilst the finite temperature discus- 
sion of Section VI is based on the single site phase shifts 
and the so-called L(ocal)-SIC implementation^ in terms 
of multiple scattering theory, briefly outlined in the fol- 
lowing section. 



III. LOCAL SELF-INTERACTION APPROACH 

In this local formulation of SIC we use the multiple 
scattering formalism and concentrate our attention on 
the phase shifts of electrons scattering from ions in a 
solid. If a phase shift is resonant it is reminiscent of a 
bound state at positive energies, i.e., above the zero of 
the potential which in this case is the muffin-tin zero. 
The energy derivative of the phase shift is related to 
the Wigner delay-time. If this is large the electron will 
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spend a long time on the site. Such 'slow' electrons will 
be much more affected by the spurious self-interaction 
and therefore should see an Si-corrected potential. Thus 
when the phase shift has a resonance we calculate the 
self-interaction correction in this (I, m) angular momen- 
tum channel. This is accomplished by calculating the 
one-electron charge density for this channel, which de- 
fines the charge density for the self-interaction correction 
potential to be added to the LSD potential. Then the 
phase shifts using the total (SIC-LSD) potential are re- 
calculated. This is applied m channel by m channel for 
a particular angular momentum I. Like in the case of 
band by band SIC implementation, the minimization of 
the total energy determines the optimum configuration 
of (l,m) channels for self-interaction correction. There- 
fore we can associate with each of the m channels two 
potential functions, V^ c ~ LSD (r) and V^| D (r). If the to- 
tal energies of these scenarios are sufficiently close, one 
can envisage dynamical effects playing an important role 
as a consequence of tunneling between these states. 

Due to the multiple scattering aspect of this approach 
we can easily calculate Green's functions and from them 
various observables for making contact with experiments. 
Another advantage is that it can be easily generalized to 
include the coherent potential approximation, 22 23 ' 24 - 25-26 
extending the range of applications to random alloys. In 
addition, one can use it to treat static correlations be- 
yond LSD by studying pseudoalloys whose constituents 
are composed of e.g. two different states of a given sys- 
tem: one delocalized, described by the LSD potential, 
and another localized, corresponding to the SIC-LSD po- 
tential. Combined with the disordered local moments 
(DLM) formalism for spin-fluctuationsj2L2£ this allows 
also for different orientations of the local moments of the 
constituents involved. In addition, in this formulation we 
can generalize the SIC approach to finite temperature, T, 
to study its effect on the electronic total energies, E tot , 
and the electronic contribution to the entropy. 

To fully take into account the finite temperature ef- 
fects, we calculate the free energy of a (pseudo) alloy, as a 
function of temperature, volume, V, and concentration, 
c, namely 

F(T,c,V) = E tot (T,c,V)-T(s cl (T,c,V) 

(c) + S mag (c) + S vib (c)). (7) 

Here S e \ is the electronic (particle-hole) entropy, S mlx 
the mixing entropy of the pseudoalloy, S m ag the magnetic 
entropy, and SVib the entropy originating from the lattice 
vibrations. 

The full extent of the L-SIC approach is ilustrated in 
Section VI where we discuss the famous Ce a — > 7 phase 
transition. 



IV. LOCALIZATION /DELOC ALIZ ATION AND 
VALENCE TRANSITIONS IN RARE EARTH 
COMPOUNDS 

A. Valencies of rare earths, rare earth nitrides and 
rare earth sulphides 

The change of the rare earth valence as a function of 
atomic number is one of the outstanding physical prop- 
erties of elemental rare earths. Among the materials that 
exhibit valence changes are the rare earth elements and 
their sulphides where changes from trivalent to divalent 
are observed in SmS, Eu, EuS, Yb and YbS, whilst no 
valence transition is observed in the rare earth nitides. 
The lattice constants of all pnictides, with the exception 
of Ce, behave continuously. In the sulphides (Fig. ^) 
changes of the lattice constants as large as 10% occur in 
the middle and at the end of the series. These are associ- 
ated with a change in valence from trivalent to divalent in 
the sulphides whilst the pnictides remain trivalent. The 
discontinuous behaviour of the CeN lattice contant is due 
to tetravalency. The calculated SIC-LSD energy differ- 
ences between divalent and trivalent configurations are 
presented in Fig. [5] One sees that at the beginning of 
the series the rare earths, their nitrides and sulphides are 
very solidly trivalent. This tendency for trivalency de- 
creases towards the middle of the rare earth series where 
the curve more or less repeats itself from Gd onwards. 
The nitrides remain trivalent throughout the rare earth 
series whereas SmS, Eu, EuS, Yb and YbS become diva- 
lent. From Fig. [3 it is clear that the rare earth nitrides 
are most trivalent and the rare earth sulphides are least 
trivalent. The elemental rare earths are intermediate be- 
tween the nitrides and sulphides. Whilst the trends are 
well reproduced, there is an overall tendency to overem- 
phasize divalency. Thus all the results n Fig. [SJhave been 
calibrated by 43 mRy to agree with the observed valence 
transition pressure of 6 kbar in SmS. 2 

Having firmly demonstrated that the behaviour of the 
lattice constants seen in Fig. ^is caused by the valence 
behaviour, we can endeavour to ask why. In Fig. [3] the 
number of /-band electrons in the trivalent states is plot- 
ted. The crossover from the trivalent to divalent state oc- 
curs when the number of /-band electrons is above 0.7. 
At this point the localization energy wins over the band 
formation energy and it becomes more advantageous to 
localize an extra / electron. For the nitrides the number 
of occupied itinerant /-electron states stays well below 
0.7 due to the strong hybridization with the nitrogen 
p-states. The nitrides remain trivalent throughout the 
rare earth series. With respect to this the application 
of the SIC-LSD is different between atoms and solids. 
In the case of atoms the levels are occupied in steps of 
one electron and the self-interaction is always applied. 
In the solid state the occupancy of a state can vary be- 
tween zero and one and whether or not to apply the self- 
interaction is determined by a competition between the 
band formation energy (no self-interaction) and the lo- 
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FIG. 1: Lattice constants of rare earth sulphides (after Ja- 
yaraman in Reference 0). 



calization (application of the self interaction) of each of 
the bands. This leads to a definition of valence since the 
self interaction corrected localized states are well sepa- 
rated from the valence states and no longer available for 
bonding and band formation. In the following we will 
discuss some further applications for the mono-pnictides 
and mono-sulphides of Ce, Sm, Eu and Yb. 



B. Ce monopnictides and monochalcogenides 

Cerium and cerium compounds attract considerable 
attention due to the intricate electronic properties re- 
lated to the Ce /-electrons. A variety of phenom- 
ena like heavy-fermion^aaaa mixed-valence 3 ^ and Kondo 
behaviourist is encountered. The cerium monopnic- 
tides and monochalcogenides have peculiar properties as 
a function of applied hydrostatic pressure (see Table Q}. 

With the SIC-LSD approach^, we have studied the 
high pressure behaviour of CeP, CeAs, CeSb, CeB i 34 i 35 , 
CeS, CeSe and CeTe 3 ^. For all these compounds we 
have found the trivalent Ce state to be the groundstate. 
The experimental lattice constants are reproduced within 
1% accuracy. The results regarding pressure transitions 
are summarized in Table HI which contains all theoret- 
ically computed phase transition pressures and relative 
volumes on the low and high pressure sides of the tran- 
sition (given as a fraction of the equilibrium volume at 
P = 0). 

The total energy as a function of volume for CeP is 
calculated for each of the phases Bl and B2 and with the 
/-electron treated as either delocalized (normal band pic- 
ture, as implemented with LSD) or localized (SIC-LSD). 
The results are shown in Fig. 0] The lowest energy is 
found in the Bl phase with localized /-electrons and with 



FIG. 2: Energy differences (in eV) between the divalent and 
trivalent configurations for the rare earths, their sulphides and 
nitrides^*- The dashed line shows the 'experimental' values 
for the rare earth metalsi- The circles, squares and triangles 
show the calculated values for the rare earth metals, the rare 
earth sulphides and the rare earth nitrides, respectively. The 
open circles and the crosses show the calculated values for the 
rare earth metals and the rare earth sulphides, respectively. 
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FIG. 3: The correlation between rare earth valence and the 
number of band-like /-electrons. Positive values on the y- 
axis mean trivalency. Negative values on the y-axis mean 
divalency. Note that for more than 0.7 /-band electrons the 
divalent configuration becomes more favourable. 



a specific volume of 348 a 3 ) per formula unit, which co- 
incides with the experimental equilibrium volume. The 
Bl phase with delocalized /-electrons has its minimum 
at a considerably lower volume, due to the significant 
/-electron band formation energy providing a large neg- 
ative component to the pressure. From the common tan- 
gent a phase transition is predicted at a pressure of 71 
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kbar with a volume collapse of AV/Vq — 8% (change in 
volume relative to the zero pressure equilibrium volume) , 
which is in excellent agreement with the transition seen 
at 55 kbar^ The B2 structure is not as favorable for the 
CeP compound, since the calculated energy is substan- 
tially higher than that of the Bl structure. This holds 
for both localized and delocalized /-electrons. From Fig. 
01 we conclude that the B2 structure with localized /- 
electrons is never reached in CeP, while at high pressure 
a second phase transition to the B2 structure with de- 
localized /-electrons is found. The transition pressure is 
calculated to be 113 kbar and the volume collapse 12%, 
while experimentally the B1-^B2 phase transition is seen 
at 150±40 kbar^S. The experimental volume collapse is 
11 %. Some uncertainty is associated with comparing the 
total energies calculated for the Bl and B2 phases in the 
present approach. Therefore we have also investigated 
the high pressure transition in CeP with the full-potential 
LMTO method^, where no problem of this sort occurs. 
In this case we have found a transition pressure of 167 
kbar, which is somewhat higher, but in good agreement 
with the experimental value. 

The results reported in Table U show that all of the 
observed pressure transitions in the cerium pnictides and 
chalcogenides are indeed reproducedi242^& The total en- 
ergy curves look rather similar to those of CeP in Fig. 0] 
but minor changes in the relative positions occur when 
the ligand is varied. The localized phases are generally 
more favored when the ligand ion becomes heavier, and 
as a consequence, in CeAs no isostructural derealization 
transition occurs in the Bl structure. Instead a tran- 
sition directly from the Bl structure with localized /- 
electrons to the B2 structure with delocalized /-electrons 
occurs, in agreement with experiment. In CeSb and CeBi 
the first high pressure transition to occur is from Bl to 
B2, with localized /-electrons in both cases, and only at 
higher pressures is a derealization transition predicted 
to take place. The calculated transition pressures are 
only slightly above the ranges studied experimentally. In 
this work only the B2 structure was considered for the 
second transition, but in reality the valence transition 
which eventually must take place in CeSb and CeBi may 
involve another high pressure phase. 

In CeS the first transition occurs to the Bl phase 
with delocalized /-electrons, i.e., the theory predicts 
an isostructural phase transition in CeS. The calculated 
transition pressure is 101 kbar with a volume collapse of 
6%. These findings are in excellent agreement with the 
experiment of Ref. 40, but at variance with the results 
of Ref. l4ll where no discontinuity in the pV-cxave is 
observed. These results may indicate the proximity of a 
critical point. At higher pressures CeS transforms into 
the B2 phase. According to the present calculations this 
occurs in two steps. First, at a pressure of 243 kbar, 
CeS goes into the trivalent B2 phase with a 4.6% volume 
change. In the second step, at a pressure of 295 kbar, 
the tetravalent B2 phase is reached with a 3.6% volume 
collapse. Thus, CeS reenters the localized regime for a 



very narrow pressure range. Some caution is necessary 
before accepting this rather peculiar behaviour. First 
of all the present calculations have been performed at 
T = K, and thermal fluctuations might easily merge 
the two high-pressure transitions into one. Secondly, the 
SIC-LSD calculational scheme only implements rather 
idealized pictures of either completely localized or com- 
pletely delocalized Ce /-electrons. Most likely, the Ce /- 
electrons enter into a complicated Kondo screened state 
in the B2 phase, which would alter the energetics of the 
B2 phase in a way we are currently not able to address. 
One could speculate that the Kondo screening at T = 
K would interpolate between the ideal localized situation 
at large volumes and the ideal delocalized situation at 
small volumes. Depending on the details of this screen- 
ing it could turn the transition into a continuous valence 
transition. Unfortunately, no experiments have been per- 
formed beyond 250 kbar,« The almost certain theoreti- 
cal prediction made here is that at high pressures the B2 
phase will be reached. However, experimental verifica- 
tion is needed. Since we have found^ that our present 
calculational approach tends to underestimate transition 
pressures for transitions from the Bl structure to the B2 
structure, the experimental transition pressure for the 
B1^B2 transition in CeS would most likely be in the 
range of 300-350 kbar, which is easily within experimen- 
tal reach. Of course, we can not rule out that other 
crystal structures may become important at high pres- 
sures. In elemental cerium as well as other rare earths, 
low-symmetry crystal phases occur when the /-electron 
derealization sets in£2£ 

In both CeSe and CeTe the only pressure transition 
observed is that from Bl to B2 with localized /-electrons 
in both phases. These are also first to occur according to 
the calculations, while valence transitions are predicted 
in the range of 400 kbar. Thus, the situation here is 
quite similar to that in CeSb and CeBi, apart from the 
tetragonal distortion in these compounds which was not 
found for CeSe and CeTe^& 



C. Sm monopnictides and monochalcogenides 

One of the most studied rare earth compounds is 
g m g50 J 5i 4 52 i 53 j QW temperature and zero pressure it 
crystallizes in the NaCl structure exhibiting a semicon- 
ducting behaviour. At a moderate pressure of ~ 6.5 kbar 
SmS reverts to a metallic phase with a significant volume 
collapse of 13.5%54, retaining however the NaCl struc- 
ture. Valence transitions can also be brought about by 
alloying the SmS lattice with other trivalent ions, such as 
Y, La, Ce or Gd^2^ Similar valence instabilities are ob- 
served in SmSe and SmTe^i^SLS, which also crystallize 
in the NaCl structure. For these compounds the volume 
changes continuously, but anomalously, with pressure (at 
room temperature From the photoemission studies 
it is concluded that SmSe and SmTe at amibient pres- 
sure, like SmS, are also of predominantly divalent / 6 
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FIG. 4: Cohesive energy of CeP (in Ry/formula unit) as a 
function of specific volume (in do /formula unit). Two crys- 
tal structures, the Bl and the B2, are considered, and each 
with two different treatments of the Ce /-electrons. The full 
drawn curves correspond to calculations with one localized /- 
electron per Ce atom, while the dashed curves correspond to 
itinerant /-electrons. The dotted line marks the common tan- 
gent at the isostructural phase transition in the Bl structure. 



character— while the monopnictides show clear signals 
of pure trivalent / 5 ions^SS. 

Figure shows the calculated lattice constants 
and valence stabilities for the Sm pnictides and 
chalcogenidesSiiS. Specifically, the total energy differ- 
ence is calculated for the scenarios of trivalent and diva- 
lent rare earth ions. A positive energy difference implies 
that the divalent configuration is preferred. For the Sm 
pnicides, the calculations reveal a strong preference for 
the f 5 configuration in the early pnictides, with the en- 
ergy difference of 1.8 eV per formula unit in SmN. For 
the heavier Sm pnictides, the / 6 configuration becomes 
more and more advantageous, and for Bi it is only 0.08 
eV higher than the trivalent configuration. This is the 
same trend towards more localised phases, as seen in the 
Ce compounds of the previous section, when the ligand 
ion becomes heavier. 

Moving to the Sm chalcogenides, already in the Sm 
monoxide the / 6 configuration is found to be most favor- 
able, by 0.08 eV, and in SmS by 0.20 eV. Hence, the SIC- 
LSD total energy predicts a valence transition of Sm be- 
tween the Sm pnictides and the Sm chalcogenides. This 
is not in complete agreement with the experimental pic- 
ture, according to which the divalent and intermediate- 
valent states are almost degenerate in SmS, while SmO 
is trivalent and metallic 63-64 . Thus, it appears that the 
SIC-LSD total energy functional overestimates the ten- 
dency to form the divalent configuration of Sm, by ap- 
proximately 15 mRy, in SmS. Assuming a similar error 
for all Sm compounds, this would imply that the calcu- 
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FIG. 5: a): Trivalent-divalent energy difference, AE = 
E(f 5 ) — E(f 6 ), of samarium compounds. Positive AE means 
divalent, negative AE means trivalent. The dashed line marks 
the calibrated curve—, b): Comparison of experimental and 
theoretical lattice constants of Sm mono-pnictides and mono- 
chacogenides compoundsi— Experimental values are marked 
with solid circles, while lattice constants calculated assum- 
ing a divalent (trivalent) Sm ion are marked with stars (dia- 
monds). 



lated energy balance curve in Fig. [S] should be lowered 
by approximately 15 mRy. The dashed line of the figure 
shows the energy difference with such a correction. This 
switches the balance in favor of trivalency for SmO, in 
accord with experiments^ 3 * 6 ^. Note that this is a reduc- 
tion, due to the inclusion of the spin-orbit interaction, 
on the 43 mRy calibration energy quoted earlier and ap- 
plied in Fig. |21 In other words, the spin-orbit interaction 
accounts for 28 mRy of the 43 mRy calibration energy. 
The remaining 15 mRy is therefore /-/ correlations not 
accounted for by the SIC-LSD. The lattice constants are 
seen to be in excellent agreement with experimental value 
for all compounds, corroborating the conclusion that a 
valence shift occurs between SmO and SmS. 

The trivalent phase of the chalcogenides becomes rele- 
vant at high pressure. In this phase the localised f 5 Sm 
ions coexist with a partly occupied narrow /-band, ef- 
fectively describing an intermediate valent phase^i. The 
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FIG. 6: The values of the lattice constant for the europium 
chalcogenides from experiment!— (triangles) and calculation— 
(crosses). Also shown is the calculated energy difference be- 
tween the divalent and trivalent form (squares) . The positive 
values mean divalent europium chalcogenides. 
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FIG. 7: The values of the lattice constant for the europium 
pnictides from the SIC-LSD calculation. The circles represent 
the calculated values^ and the triangles are the experimental 
values p— Also shown is the calculated energy difference be- 
tween the divalent and trivalent form. Negative values mean 
trivalent, positive values mean divalent. 



calculated and measured transition pressures are listed in 
Table HH The good agreement both for transition pres- 
sures and volume collapses proves that the bonding of the 
high pressure phase is well described in the SIC-LSD ap- 
proximation, even if the true many-body wavefunction of 
the intermediate valence phase is much more complicated 
than the corresponding SIC-LSD wavefunction. This is in 
line with the general philosophy of the density functional 
approach of obtaining good total energy estimates from 
simple reference systems (non-interacting electrons) . The 
present theory cannot describe the continuous nature of 
the transition observed for SmSe and SmTc. The exper- 
iments were all conducted at room temperature and it 
would be interesting to investigate whether the continu- 
ous transition would exist at low temperature as well. 



D. Eu monopnictides and monochalcogenides 

Europium chalcogenides and most of the pnictides 
crystallise in the simple NaCl crystal structure and hence 
form a series that can be studied within first principles 
theory relatively easily. Recently the chalcogenides have 
attracted a lot of attention due to their potential appli- 
cations in spintronic and spin filtering devices— i 

The SIC-LSD method is applied to study the electronic 
structure of the Eu compounds in the assumed ferromag- 
netic state in both the divalent and trivalent configura- 
tions. Fig. shows the calculated and experimental lat- 
tice constants for all the chalcogenides r 6 - Also shown are 
the energy differences between the two valence states. It 
is clear that all the europium chalcogenides are divalent. 



This is as expected from simple shell- filling grounds. It 
can also be seen from this figure that the calculated lat- 
tice constants are in good agreement with experiment. 
The energy difference between the two valence states is 
fairly independent of chalcogenide for S, Se and Te and 
one can observe that there is a clear correlation between 
the lattice constant and the difference in energy between 
the divalent and trivalent states. 

The calculated transition pressures for EuO and EuS 
are compared with experiment in Table IIIII The lower 
transition pressure and smaller volume in EuS compared 
to EuO is reproduced by the theory. Experimentally, 
the transition of EuO (at room temperature) is continu- 
ous, which the present theory cannot describe. For EuS 
the experiments show no anomalous compression curved, 
but the band gap closes at 160 kbar, just before the struc- 
tural transition to the CsCl structure (at 200 kbar)2£. 
However due to the LMTO-ASA approximations, signifi- 
cant uncertainty persists in the values of the total energy 
differences between different crystal structures. Also the 
spin-orbit interaction can significantly alter the results: 
we found that without spin-orbit the structural transition 
occurs at 137 kba»iS, however without an isostructural 
transition occurring first. 

Fig. [3 shows the equilibrium lattice constants for the 
europium pnictides. Also shown are the energy differ- 
ences between the divalent and trivalent configurations. 
It is clear that EuN, EuP and EuAs are trivalent and 
EuSb is divalent. These results are in full agreement with 
Hulliger 72 who states that Eu ions in the EuN and EuP 
are known to be trivalent, while EuAs is known to contain 
some divalent ions. We are not aware of any definitive 
measurement of the valence state of EuSb, but clearly 
if the trend continues it will be divalent. Although the 
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calculation predicts that EuAs, with the rocksalt crys- 
tal structure, is trivalent, EuAs has in reality the Na2C>2 
crystal structure. This is a distortion of the NiAs struc- 
ture due to the formation of anion-anion pairsi 2 . We 
speculate that the occurence of Na2C>2 crystal structure 
and divalent ions is intimately connected. These diva- 
lent ions do not occur in the rocksalt crystal structure 
and hence explain the theoretical underestimate in the 
value of the lattice constants as seen in Fig. The pres- 
ence of divalent ions in EuAs would also make the valence 
energy difference curve of Fig. 0more continuous. These 
results are qualitatively similar to those obtained for the 
ytterbium pnictides2i2^ (see below) where there is also 
an increasing tendency for divalency as we go down the 
pnictidc column of the Periodic Table. However, in that 
case the divalent state is not reached. 



E. Yb monopnictides and monochalcogenides 
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Here the application of the SIC-LSD method to Yb and 
a number of its compounds is discussed J 4 ' 7 ^ We concen- 
trate in particular on the valence of Yb ion in these sys- 
tems, and some of the valence transitions, in particular 
in YbS. Fig. [5] shows for Yb compounds the / electron 
difference between the divalent and trivalent configura- 
tions versus their energy difference. The behaviour is 
close to linear. The closest we come to change of one elec- 
tron occurs in YbN, specifically a change of 0.8 electrons, 
for which also the highest energy difference between the 
two electronic configurations is obtained. This maximum 
electron difference of 0.8 gradually decreases to approx- 
imately 0.5 for YbPd, YbSb, YbBi and YbBiPt and for 
the four divalent Yb systems we obtain between 0.22 to 
0.35 electrons. For these divalent systems the difference 
in / electron count has nearly disappeared and this seems 
to herald the arrival of the divalent behaviour. 

Our calculations can also make contact with the 
pressure-volume measurements which e.g. for YbS indi- 
cate anomalous behaviour around f 00 kbar which could 
be associated with intermediate valence 7 ^. In particular, 
for this system a trivalent state cannot seemingly be re- 
alized. In this case, our study could shed some light on 
properties of this intermediate valence state. From the 
common tangent construction of the total energies as a 
function of volume for both divalent and trivalent YbS, 
we obtain with the SIC approach a transition pressure 
which is ~75 kbar and agrees well with the experimental 
value of about ~1 00 kbar (Fig. EJ). As seen in Fig. the 
experimental anomaly is not really correlated with an in- 
teger change in valence, but is as a matter of fact due to 
an / electron derealization. The quantitative theoretical 
description of this seemingly continuous valence transi- 
tion calls for a more elaborate theory than presented here. 
The small change in / electron occupancy of 0.3 electrons 
found with the SIC approach upon derealization of an 
/ electron suggests the occurence of intermediate valence 
of 2.3. This is consistent with the experimental estimate 
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FIG. 8: Divalent-trivalent energy difference (in eV) of Yb 
compounds versus the total /-occupancy difference in the di- 
valent and trivalent configurations, as given by the SIC-LSD 
approach. 



of an intermediate valence of 2.4i2i 



V. VALENCE TRANSITIONS IN ACTINIDE 
COMPOUNDS 

A. Elemental actindides 

Compared to the rare earth 4/-states, the 5 /-states in 
the actinides are less inert, and in Th, Pa, and U, play an 
active role in the cohesion, as manifested by the low sym- 
metry crystal structures, low specific volumes, and large 
bulk moduli. A localization transition occurs when going 
from Pu to Ampi and in the later actinides, from Am 
to Es, the /-electrons are non-bonding, high-symmetry 
crystal structures are attained, the specific volumes are 
large and the bulk modulii relatively small. Pu is situated 
at the borderline between these two competing pictures, 
and its very complex phase diagram implies that the /- 
electron properties are of particularly intricate nature. 
The low-symmetry cv-Pu ground state is well reproduced 
by LDA calulations, 78 whilst the large-volume fee phase 
of <5-Pu, is believed to be characterized by localized /- 
electrons. 

The SIC-LSD approach was applied to the calculation 
of the total energies as a function of atomic volumes for 
the actinide elements from U to Fm. 10,79 Ferromagnetic 
and paramagnetic arrangements were considered in the 
fee structure, while antiferromagnetism was investigated 
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FIG. 9: Equation of state for YbS as calculated by the SIC- 
LSD method. The two theoretical curves correspond to 14 
(solid line) and 13 (dashed line) localized / electrons, respec- 
tively, while the dots are the experimental data of Ref. l76l 
The dotted line marks the theoretical transition. 



in the hep structure. The overall total energy minima 
were found to occur for the trivalent configuration in the 
case of Am to Cf, while Es and Fm have their minima 
for the divalent state. Experimentally, Es is believed to 
be divalent, due to its large lattice constant, while Fm 
has never been prepared in the solid state. The U and 
Np metals have the lowest energies when the /-electrons 
are fully delocalized. The position of Pu at the crossover 
between localized and delocalized /-electron behaviour is 
revealed in the very different outcome of the calculations 
depending on the magnetism. In a magnetic calculation, 
the trivalent configuration has the lowest energy, how- 
ever with a volume ~ 30% higher than the experimental 
volume of £-Pu, while a paramagnetic treatment leads to 
virtual degeneracy of any of the localization scenarios /° 
to / 4 . The trends of the energy difference between the 
divalent and trivalent configurations are qualitatively the 
same in the actinides and the lanthanides and are gov- 
erned by the relative stability of the half-filled / 7 shell. 
In Table IIVI the calculated equilibrium volumes, bulk 
modulii, and transition pressures for onset of derealiza- 
tion are shown. Overall, the experimental data are well 
reproduced by the SIC-LSD calculations. 

The peculiar role played by Pu in the series of elemen- 
tal actinides is illustrated in figure which shows the 
total energy as a function of volume for all localization 
scenarios from f° to / 6 in a paramagnetic treatment. Re- 
markably, within 0.03 eV/atom the scenarios with 0,1,2,3 
or 4 localized /-electrons are degenerate, while localizing 
5 or 6 /-electrons is less favourable. The interpretation 
of this is that Pu at zero pressure exists in a complex 
quantum state vividly fluctuating between localized and 
delocalized /-electrons. This is in accord with the ex- 
perimental observation of several allotropic forms of Pu, 



including the low-volume a-phase and the high-volume 
(5-phase, the latter being only stable at elevated temper- 
ature or by alloying. In Table ITvl we compare the data of 
the / 4 state to the 8 phase. If magnetic order is imposed, 
in either ferro- or antiferro-magnetic arrangement, the / 5 
configuration becomes the ground state, i.e. formation of 
an ordered magnetic moment favours localization. The 
too large a volume associated with this solution indicates 
that this localization tendency is in fact overestimated, 
as the SIC-LSD formalism does not account properly for 
the strong quantum fluctuations taking place for a Pu ion 
in the metallic phase. Experimentally, there is no sign of 
magnetic moments in elemental Pu. 83 



In Fig. ^2 the total energy as a function of volume is 
shown for several /" configurations in Am. As already 
mentioned, the trivalent / 6 configuration has the low- 
est energy, with an equilibrium volume of 201 a.u., while 
the experimental volume is 198 a.u.. Upon compression 
the band formation energy increases and the energy dif- 
ference between / 5 and f 6 decreases. At a volume of 
about 140 a.u. the / 5 and / 6 total energy curves cross, 
and from then on the /-electron band states start to 
dominate. Around V — 100 a.u. the /-derealization 
is complete. This is in good agreement with the ex- 
perimental observation of a structural phase transition 
(Amll-Amlll) taking place at 100 kbar and at a compres- 
sion corresponding to 77% of the equilibrium volume^ 
which is interpreted as the onset of /-bonding in Am. 
A rough theoretical estimate of the transition pressure, 
given by the slope of the E(V) curves at the crossing 
point, leads to a value of 160 kbar. This estimate con- 
stitutes an upper bound to the derealization pressure, 
since the high-valence SIC-LSD configurations represent 
a rather limited many-body wavefunction for a correlated 
state with actively bonding /—electrons. In addition, the 
correct (orthorhombic) high pressure phase was approx- 
imated here by the fee structure, and finite temperature 
may also lower the phase transition pressure. Experimen- 
tally, the occurrence of the AmIV phase at V/Vo = 0.63 
at p = 175 kbar probably marks the endpoint of the /- 
electron derealization process in Am. 



Figurc lT21 shows the upper bounds for a transition pres- 
sure to a phase with bonding /—electrons. The agree- 
ment with experiment varies accross the series, but is 
generally not bad (within a factor of 2). This qualita- 
tive agreement shows that the magnitude of the SIC- 
LSD localization energy is roughly correct. FigurelT2lalso 
shows the comparison of the volume ranges over which 
/-electron derealization occurs for the actinide elements 
Pu, Am, Cm and Bk, as given by this theory and ex- 
periment. The qualitative trend as given by the relative 
stability of Cm(/ 7 ), is well reproduced by the calcula- 
tions, with a tendency of the calculated end-points of 
/-transition to be at too low a volume. 
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FIG. 10: Total energy of Pu in the paramagnetic fee state. 
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FIG. 11: Total energy versus volume for Americium in the 
paramagnetic fee structure. Several configurations of the lo- 
calized 5/ subshell are considered. 



B. Actinide monopnictides and monochalcogenides 

In the actinide monopnictides and monochalcogenides, 
which all crystallize in the NaCl structure at ambient 
conditions, the actinide-actinide separations are larger 
than in the elemental metals, and the tendency towards 
/-electron localization can already be observed from Np 
compounds onwards. Figure El displays the calculated 
SIC-LSD ground state configurations through the series 
of U-, Np-, Pu-, Am-, and Cm-mono-pnictides and - 
mono-chalcogenides. The calculations reveal clear trends 
towards more and more actively bonding /-electrons for 
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FIG. 12: Trends in the ranges of volumes of / electron de- 
localization in the actinide elements. Onset (completion) of 
derealization is marked with balls (squares), for each element 
with experimental data (blue and full line) to the left and the- 
oretical data (red and dashed line) to the right. Volumes are 
given relative to the (localized) equilibrium volume at ambient 
conditions. Experimental ranges are defined by the smallest 
volumes observed in the high-symmetry (fee) phase and the 
largest volumes observed for the low-sym metry (q-U type) 
phases in high pressure experiments (Refs. |82|84||85l for Pu 
the range is defined by the zero pressure volumes of the 3- 
and a-phases). The theoretical ranges are calculated within 
the fee structure only. 



a) lighter actinides, and b) lighter ligands. For the lighter 
actinides, the /-orbitals are more extended, leading to 
larger overlaps with their nearest neighbours, and smaller 
self-interaction corrections, both of these effects favour- 
ing band formation. For the lighter ligands, in particular 
N and O, both the volume is decreased and ionicity is 
larger, the first of these effects leading to larger direct 
actinide-actinide overlap, and the latter effect favouring 
charge transfer. 

The Cm compounds are the most localized systems, 
all exhibiting Cm in the trivalent f 7 configuration. The 
f 7 shell is so stable that variations of the ligand can- 
not disrupt its stability, and scenarios with either one 
more or less localized /-electron have distinctly higher 
energies. Trivalency prevails in the Am-compounds, but 
the stability of the f 7 shell causes the divalent Am state 
to be important in AmTe and AmPo. In the Pu com- 
pounds, the trivalent state also dominates, but for the 
lighter ligands /-electron derealization sets in. In the 
Np compounds the tetravalent state dominates, while in 
the U compounds pentavalent states occur for the lighter 
ligands. 

The density of states (DOS) of the actinide arsenides 
are shown in Fig. 1141 with both trivalent and tetravalent 
actinide ions. In the trivalent case, the non-(SIC) /- 
degrees of freedom give rise to narrow unoccupied bands 
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FIG. 13: Trends in localization through the AcX series. For 
each actinide, Ac=U, Np, Pu, Am, Cm, a block of 10 lig- 
ands are considered; the pnictides X=N, P, As, Sb, Bi, and 
the chalcogenides X= O, S, Se, Te, and Po. The numbers 
designate the calculated Ac valence for that particular AcX 
compound. Where two numbers are given, the corresponding 
valences are degenerate. 



5/ electrons^! and in UPd3 this same shoulder is seen 
about 1 eV below the Fermi level, and is interpreted in 
terms of localized / electrons^ Sharp features, usually 
associated with localized /-electrons, are not observed. 
In UPt3_ K Pda; alloys, gradually substituting Pt by Pd, 
leads to a transition from the pentavalent to tetrava- 
lent groundstate between UPt2Pd and UPtPd2, as can 
be seen from Figs. 115b and 115b respectively. An addi- 
tional /-electron gets localized, a transition that appears 
to be driven mainly by the changes in the electronic struc- 
ture related to the decreasing fd hybridization when re- 
placing Pt by Pd. The opposite effect is obtained by 
putting UPd3 under pressure. The total energy calcula- 
tions for UPd3 show that the pentavalent configuration 
becomes energetically favourable at pressures of approx- 
imately 250 kbar. From the observed similarities in the 
pentavalent DOS of both UPd3 and UPt^^i, and given 
that UPt3 is a heavy fermion material at zero pressure, 
we expect UPd3 to become heavy fermion under pressure. 
Experimental efforts 90 to see the proposed transtion of 
UPd3 under pressure failed. This can be either due to 
an underestimate in the present theory of the transition 
pressure, or finite temperature effects in the experiment. 



above the Fermi level. In the tetravalent case the addi- 
tional delocalized /-electron appears as an extra /-band. 
In Cm, this band appears far below the Fermi level, while 
in Am, Pu and Np this band lies just below the Fermi 
level. The band formation due to this extra band is suffi- 
ciently large in NpAs to outweigh the localization energy, 
and the tetravalent configuration becomes the ground 
state. 



C. Uranium compounds UPt3 and UPd3 

Using the SIC-LSD method, we have similarly studied 
the 5 / valence configurations of U in UPt3 and UPd3i2& 
Both compounds are isoelectronic as far as the valence 
electrons are concerned, however the 5d electrons are 
less tightly bound to the nucleus than the Ad electrons, 
which results in a increased overlap of the U / and the 
transition metal d orbitals in UPt3. Correspondingly, 
the total energy calculations determine the tetravalent 
/ 2 and the pentavalent f 1 configurations for UPd3 (Fig. 
ITok ) and UPt3 (Fig. 1X31 1) respectively. The calculated 
equilibrium volumes, VuPda ~ 474.1(a.u. 3 ) and Vupt 3 = 
469.7(a.u. 3 ), compare well with the experimental values 
of V^ d = 469.5(a.u. 3 ) and Vg% a = 472.9(a.u. 3 )S In 
both UFM3 and UPt3, the small energy differences be- 
tween the tetravalent and pentavalent configurations in- 
dicate near degeneracy (within ^5 mRy), rather than 
a clearly preferred / 2 groundstate for UPd3 and an f 1 
groundstate for UPt3. This in a way mirrors the results 
from XPS measurements, where in UPt3 a weak shoul- 
der at the Fermi energy indicates the presence of itinerant 



VI. VALENCE TRANSITION IN CE AT FINITE 
TEMPERATURE 

Cerium is the first element in the Periodic Table that 
contains an / electron, and shows an interesting phase 
diagram. 91 In particular, the isostructural (fee — * fee) 
a — 7 phase transition is associated with a 15%-17% 
volume collapse and total quenching of the magnetic 
moment & The low-pressure 7-phase shows a local mag- 
netic moment, and is associated with a trivalent config- 
uration of Ce ion. At the temperatures in which the 
7-phase is accessible, it is in a paramagnetic disordered 
local moment state. Increasing the pressure, the material 
first transforms into the a-phase, which is indicated to 
be in an intermediate valence state with quenched mag- 
netic moment. At high pressures (50 kbar at room tem- 
perature) Ce eventually transforms into the tetravalent 
a'-phase. With increasing temperature, the a — 7 phase 
transition shifts to higher pressures, ending in a critical 
point (600K, 20 kbar), above which there is a continuous 
crossover between the two phases. 

To describe the full phase diagram of the Ce a— 7 phase 
transition we have modelled Ce as a pseudoalloyjS 2 ^ in 
the spirit of the Hubbard III approximation^, consisting 
of the trivalent (SIC-LSD) Ce atoms with concentration 
c, and the tetravalent (LDA) Ce atoms with the concen- 
tration (1-c). In order to properly take into account the 
disordered local moments of the trivalent Ce atoms in the 
7-phase, equal probabilities for up and down orientations 
of their local moments were adopted. Assuming homo- 
geneous randomness, this ternary pseudoalloy can be de- 
scribed by the coherent potential approximation (CPA). 
The respective concentrations of the trivalent and the 
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FIG. 14: Densities of states (in states/(Ry formula unit)) for the actinide arsenides NpAs, PuAs, AmAs, CmAs, with a trivalent 
(left), or tetravalent (right) actinide ion, respectively. The solid and dotted lines represent the / projected- and total densities 
of states, respectively. The energies are given in Ry, with the Fermi level at energy zero. 



tetravalent Ce ions in the pseudoalloy are then deter- 
mined by minimizing the free energy for each volume 
and temperature with respect to the concentration c. 

The free energy of the physical system at a given vol- 
ume can be obtained by evaluating the concentration 
dependent free energy at the minimizing concentration 

Cmin ■ 

F(T, V) = F(T, Cmin, V) . (8) 

These free energies are displayed in Fig. El which clearly 
shows a double- well behaviour, corresponding to the two 
separate phases of Ce, at low temperatures, which is 
gradually smoothened out with increasing temperatures. 
Furthermore one finds that, at elevated temperatures, 
the free energy is mainly lowered for large lattice con- 
stants, corresponding to the 7-phase, with its larger en- 
tropy. Inserting the minimizing concentration c m ; n into 
the pressure- volume relation 

p(T, V) = p(T, c min , V) = —^F(T, c min , V), (9) 

allows to calculate the isotherms of Ce, which are dis- 
played in Fig. El It can be seen that the average valence, 
close to the coexistence line, gradually changes with in- 
creasing temperature. Above the critical temperature, 
the valence changes continuously with increasing pres- 
sure from trivalent to tetravalent Ce ion. 

In Fig. El we present the phase diagram, obtained 
from the free energies of the a-7 pseudoalloy, with the 



7-phase described by the DLM approach. It can clearly 
be seen in the figure how the transition becomes contin- 
uous above the critical temperature. The experimentally 
observed critical point (600K, 20 kbar) falls on top of 
the calculated phase separation line, which starts at the 
zero temperature transition pressure of -7.4 kbar. This 
means that the slope of the phase separation line is in 
very good agreement with experiments. The calculated 
critical temperature overestimates the experimental one 
by roughly a factor of 2, which is still reasonable con- 
sidering that the critical temperature is very sensitive to 
various small details of the calculations and in particular 
the theoretical lattice parameters of both the Ce phases. 
Note that our calculated value of 169 K for the T c at 
zero pressure compares well with the experimental value 
ofl41±10K. 

Finally we examine in more detail the discontinuity 
across the phase separation line. Fig. Elshows the mag- 
nitude of the discontinuities for the various contributions 
of the Gibbs free energy. As expected, all contributions 
vanish at the critical temperature, above which there is 
a continuous cross-over between the a- and the 7-phase. 
It also can be seen from this figure that the entropy dis- 
continuity is by far the largest contribution. The phase 
transition is therefore driven by entropy, rather than by 
energetics. The entropy discontinuity itself is mainly due 
to the magnetic entropy. Thus it is the entropy, and not 
the internal energy that drives the Ce a — 7 phase tran- 
sition. Analysis of experimental data leads to the same 
conclusion. 95 
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FIG. 15: Total Energy (in Ry/f.u.) versus volume (in (a.u.) 3 ) 
for a) UPd 3 , b) UPtPd 2 , c) UPt 2 Pd, and d) UPt 3 . 



VII. SUMMARY AND CONCLUSION 

In this article a review has been given of some ap- 
plications of the SIC-LSD to the calculation of the va- 
lences of /-electron systems. A methodology has been 
presented which is able to determine valence changes as 
a function of pressure and chemical composition. An 
important finding of this work is the discovery of the 
dual character of / electrons. The number of localized 
/-electrons defines the valence. Furthermore, a corre- 
lation between a change in valence and the number of 
band-like /-electrons has been established. 

The versatility and functionality of the SIC-LSD has 
increased with the L-SIC. In particular as the example 
of the p-T plot in Fig. 1181 shows a finite T generalization 
of the method has been successfully developed. This will 



allow us to perform the finite T study of the p-V curves 
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FIG. 16: The free energies as function of the volume for 
the temperatures (highest curve), 200, 400, 600, 800, 1000, 
1200, 1400 and 1600 K (lowest curve). The zero of energy is 
arbitrary. 

of several of the compounds discussed in this paper and 
to determine if the continuous transitions as some times 
seen experimentally can be obtained. 
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FIG. 17: Calculated isotherms for the temperatures T=0 
(lowest curve), 200, 400, 600, 800, 1000, 1200, 1400, and, 
1600K (highest curve). The color indicates the fraction of 
localized electrons: blue is all localized (7-phase) and red is 
all delocalized (a-phase). 
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FIG. 18: Phase diagram obtained for the pseudoalloy, com- 
posed of a- and 7-Ce. The crosses indicate the calculated and 
experimental critical points. 
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TABLE I: Calculated transition pressures for the electronic and structural phase transitions in the cerium monopnictides and 
monochalcogenides. Also quoted are the specific volumes (relatively to the zero pressure equilibrium volume) on the two sides 
of the transition Ref. I49I The notation (d) and (1) refers to calculations with delocalized or localized Ce /-electrons, i.e. 
tetravalent or trivalent Ce atoms. B2* denotes the distorted B2 structure. 

The subscripts refer to the following references: a: Ref. EH; b: Ref. EH c: Ref. H d: Ref. |H e: Ref. |H f: Ref. |H g: Ref. 
E3; h: Ref. |H i: Ref. |U j: Ref. H 



compound 


transition 




Pt (kbar) 


Vi/Vo 


V-t/Vo 








theo. 


expt. 


theo. 


expt. 


theo expt. 


CeN 


Bl(d) - 


B2(d) 


620 




0.760 


- 


0.724 


CeP 


Bl(l) -> 


Bl(d) 


71 


90 a ,55 b 


0.933 


0.89 a 


0.853 0.85 a 


CeP 


Bl(d) - 


B2(d) 


113 


150(40) a ,250 J 0.827 


0.82 a 


0.706 0.71" 


CeAs 


Bl(l) - 


B2(d) 


114 


140(20) C ,210 J 0.893 


0.84 c 


0.713 0.73 c 


CeSb 


Bl(l) - 


B2*(l) 


70 


85(25) d ,150 j 


0.922 


0.90 d 


0.813 0.80 d 


CeSb 


B2*(l) -> 


B2*(d) 


252 




0.717 




0.680 - 


CeBi 


Bl(l) - 


B2*(l) 


88 


90(40) e 


0.897 


0.87 e 


0.789 0.78 e 


CeBi 


B2*(l) -> 


B2*(d) 


370 




0.666 




0.638 - 


CeS 


Bl(l) - 


Bl(d) 


101 


- / ,125(15) 9 


0.918 


0.93 9 


0.855 0.88 s 


CeS 


Bl(d) - 


> B2(l) 


243 




0.788 




0.742 - 


CeS 


B2(l) - 


B2(d) 


295 




0.724 




0.688 - 


CeSe 


BIO) - 


B2(l) 


124 


170(30)' 1 


0.890 


0.86 h 


0.779 0.77 h 


CeSe 


B2(l) - 


B2(d) 


377 




0.683 




0.652 - 


CeTe 


Bl(l) - 


B2(l) 


74 


55(25)* 


0.915 


0.93* 


0.798 0.84* 


CeTe 


B2(l) -» 


B2(d) 


435 




0.647 




0.623 - 



TABLE II: Calculated isostructural transition pressures, Pt (in GPa; 1 GPa=10 kbar), and volume changes (in %), of Sm 
monochalcogenides. Experimentally, the transition of SmS is discontinuous, while those of SmSe and SmTe (at room tempera- 
ture) are continuous. 

The subscripts refer to the following references: d : Ref. I5"3 ; e : Insulator-metal transition of Ref. IB"?! Present author's 
estimates from figures of Ref. [H^ and 9 : Ref. IdH The volume changes for SmSe and SmTe are obtained by extrapolation over 
the transition range. 



Compound 




Pt (GPa) 


Volume collapse (%) 




Theory 


Expt. 


Theory Expt. 


SmS 


0.1 


0.65 d , 1.24 e 


11.1 13. 5 d , 13.8 e 


SmSe 


3.3 


~ 4 d , 3.4 e , 3 - 9 f ,2.6 - 4 9 


9.8 8 d , ll f , 7 9 


SmTe 


6.2 2 


-8 d , 5.2 e , 6-8 / , 4.6 - 7.5 9 


8.4 9 f ,7 9 



TABLE III: Calculated isostructural transition pressures, Pt (in GPa; 1 GPa=10 kbar), and volume changes (in %), of Eu 
monochalcogenides. 

The subscripts refer to the following references: a : Ref. liist b : Ref. lo9t c : Insulator-metal transition, Ref. [73 



Compound 


Pt (GPa) 
Theory Expt. 


Volume collapse (%) 
Theory Expt. 


EuO 
EuS 


19.3 SO^.IS-SO 6 
11.6 16 c 


6.3 5 a 
5.7 C 
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TABLE IV: Calculated— and experimental equilbrium volumes, V, bulk modulii, B, and /—electron derealization pressure, 
P, for the actinide elements. The magnetic ground state was used for Cm onwards, while the nonmagnetic was used for Pu 
and Am. Experimental values are from Ref. il except a : Ref. US and 6 : Ref. 





Vteo (a.u.) 


Vexp (a.u.) 


B teo (GPa) 


B eX p (GPa) 


Pt.o (GPa) 


Pexp (GPa) 


Q-PU 


123.0(-9%) 


135" 










<5-Pu 


163.8(-2.5%) 


168° 


46.0 


32 b 





~ 


Am 


200.6(+1.3%) 


198 


46.0 


45 


16 


10 


Cm 


203.8(+0.9%) 


202 


42.4 


33(5) 


70 


43 


Bk 


197.6(+6%) 


189 


35.0 


25(5) 


15 


25 


Cf 


201.8(+9%) 


185 


37.3 


49(5) 


30 


41 


Es 


256.0(-4.1%) 


267, 321 


19.5 




11 




Fm 


247.4 




29.6 




28 





